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We study the two-band effects on ultracold fermionic atoms in optical lattices by means of dynam- 
ical mean-field theory. We find that at half-filling the atomic-density-wave (ADW) state emerges 
owing to the two-band effects in the attractive interaction region, while the antiferromagnetic state 
appears in the repulsive interaction region. As the orbital splitting is increased, the quantum phase 
transitions from the ADW state to the superfluid state and from the antiferromagnetic state to 
the metallic state occur in respective regions. Systematically changing the orbital splitting and the 
interaction, we obtain the phase diagram at half-filling. The results are discussed using the effective 
boson model derived for the strong attractive interaction. 



PACS numbers: 03.75.Lm, 05.30.Fk, 73.43.Nq 



I. INTRODUCTION 

Ultracold fermionic atoms in optical lattices have at- 
tracted considerable attention. By loading fermionic 
atoms into optical lattices, diverse interaction config- 
urations can be introduced. Furthermore, Feshbach 
resonances provide the means for controlling both the 
strength of the interaction between fermionic atoms and 
its sign p| . Fascinating quantum many-body phenomena 
have been revealed by the combination of these experi- 
mental techniques. The topological change in the Fermi 
surface of K fermionic atoms was observed by increas- 
ing the band- filling 2]. By controlling the atomic inter- 
action, a band insulator in the lowest band was produced 
and then the partially populated higher bands were ob- 
served. The Mott insulating state for 40 K was realized 
quite recently by adequately tuningthe ratio between 
the interaction and kinetic energy [3|. Three features 
were identified; a suppression of doubly occupied lattice, 
a reduction of the compressibility, and the gapped mode 
in the excitation spectrum. It was argued that the re- 
sults pave the way for future studies of spin ordering. 
For 6 Li fermionic atoms with attractive interaction, a 
superfluidity of fermionic atom pairs was observed [J]. 
By increasing the depth of the lattice potential near the 
Feshbach resonance, a superfluid-insulator transition was 
observed. In the experiments it was argued that the usual 
single-band model was no longer applicable, because the 
strength of the on-site interaction exceeded the gap be- 
tween the lowest and the next-lowest bands. Accordingly, 
the effects of the higher bands have to be taken into ac- 
count for detailed investigations. 



dered states LL^I, LLJ, UM tor termiomc atoms m two- or 
three-dimensional optical lattices have been investigated 
by various theoretical methods. Recently, we have in- 
vestigated the superfluid-insulator transition of fermionic 
atoms in optical lattices beyond the conventional mean- 
field approximation, taking two-band effects into account 
(To| . It was shown that the Mott insulating state appears 
for the adequately strong interaction region at half-filling. 
For the Mott insulator in the repulsive interaction the 



fermionic atoms in each site occupy both orbitals, while 
in the attractive interaction the bosonic fermion pairs oc- 
cupy either of the two orbitals in each site. The results 
suggest that other ordered states may be possibly in- 
duced: the atomic-density-wave (ADW) state for the at- 
tractive region and the antiferromagnetic (AF) state for 
the repulsive region. However, an issue whether these or- 
dered states compete or coexist with the superfluid state 
or Mott insulating states has not yet been well investi- 
gated for two-band lattice fermionic systems. 

In this paper, we investigate the ordered state of ul- 
tracold fermionic atoms in three-dimensional optical lat- 
tices, taking the two-band effects into account. For this 
purpose, we make use of a dynamical mean-field theory 
(DMFT) [l6| , which enables us to treat local correlation 
effects precisely. We show that for the attractive inter- 
action region the ADW state is stable for a small orbital 
splitting region and that the transition to the superfluid 
state takes place with increasing the orbital splitting. For 
the repulsive interaction region, the AF state appears in 
a small orbital splitting region. As the orbital splitting 
increases, the AF state changes to the metallic state. For 
large orbital splittings, the band insulating state comes 
into existence in both regions. 

The paper is organized as follows. In Sec. II, we in- 
troduce the model Hamiltonian and explain the two-site 
DMFT method to investigate the ADW and AF states. 
In Sec. Ill, the numerical results for the attractive and 
repulsive regions at half-filling are shown. We discuss 
the most stable state for given parameters by comparing 
their energies. The results are summarized in the phase 
diagram. In Sec. IV, we discuss the numerical results 
in comparison with the effective boson model, which is 
derived for the strong attractive region. It is shown that 
both results are consistent with each other. A brief sum- 
mary is given in Sec. V. 



II. MODEL AND METHOD 



Let us consider the fermionic atoms in a optical lat- 
tice potential: V(r) = Vb(sin 2 fca; + sin 2 ky + sin 2 fcz). 
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In the low-tunneling Vq 3> E r , where E r = H 2 k 2 /2m is 
the recoil energy, each lattice potential is regarded as a 
harmonic one [T7L IHI ] . We investigate the effects of the 
lowest and next-lowest orbitals, so that the three-fold de- 
generacy of the next-lowest orbitals is neglected for sim- 
plicity. The hopping integrals between the lowest orbitals 
(t\) and between the next-lowest orbitals (£2) satisfy the 
relation t 2 ~ \/Vo/E r t%. Since Vo/E r < 10 in the exper- 
iments 0,0] 1 we approximately set that t\ = £2 =t. The 
following interactions are considered: the on-site intraor- 
bital interactions for the lowest orbital (Ui) and for the 
next-lowest orbital (U2), the interorbital interaction (W), 
and the interaction corresponding to the Hund coupling 
(J). The coupling constants of these four interactions 
satisfy the relations U 2 = (3/4) Ui and U' = J = U x /2. 
We set that U\ = U2 = U approximately and thus 
U' = J= U/2. 

The system is assumed to involve the same number of 
fcrmionic atoms in two different hypcrfine states, which 
are described as the pseudospins. The simplified model 
Hamiltonian thus obtained reads 

(i,j)aa io 

+ Uy^i^Tlial +} j (U'- J8 ata i)n i i c n i 2a', (1) 
iot iaa' 

where Ci a(T is the fermionic annihilation operator for the 



g q M = 

where p(z) is the density of states (DOS). '^a{b'),cX u> ') 
is the self-energy of the orbital a for the A(B) sub- 
lattice, which can be obtained by solving two effec- 
tive impurity models. We use a semicircular DOS, 
p{z) = A/{ttW)^1 - 4(z/W) 2 , where W is the band 
width. Since the hopping integral is assumed to be in- 
dependent of a, W — At and the DOS are the same 
for both bands. The chemical potential is set to be 
/i = U/2 + U' — J/2 so that particle- hole symmetry can 
be satisfied. In this case, two bands by the orbitals a = 1 
and 2 together are half-filling. In the following, the hop- 
ping integral t is used in units of energy. 



III. NUMERICAL RESULTS 

For the attractive interaction region U < 0, three types 
of states are considered: the superfluid state, the Mott in- 



state with pseudospin <r(=t ,1) on orbital a(—l, 2) in the 
ith lattice site, n iaa = c' iacr Ci au , and the subscript 
indicates the sum of the nearest-neighbor sites, fj, is the 
chemical potential and D is the splitting between the 
two orbitals. We assume that the intraorbital attractive 
interaction induces an s-wave superfluid state. 

For the attractive region (U < 0) we examine the 
superfluid-insulator transition, turning our attention to 
the stability of the ADW state. For the repulsive region 
(U > 0) the phase transition between the metallic state 
and the insulating state is examined, laying stress on the 
stability of the AF state. 

In DMFT, the lattice model is mapped onto a single 
impurity model connected dynamically to a heat bath. 
The Green's function is obtained via the self-consistent 
solution of this impurity problem. This retains nontrivial 
local quantum fluctuations missing in conventional mean- 
field theories. We apply here the two-site DMFT method 
[l9| , which allows us to study the Mott transitions of or- 
bitally degenerate lattice fermions qualitatively [2(1 Hl| . 
To study the superfluid of lattice fermions, we extend 
this method to the case when the superfluid order exists 

In order to examine the ADW and AF states, we divide 
the bipartite lattice into two sublattices [lfl [22j . In this 
procedure, the local Green's function has the following 
form: 

(2) 
(3) 

I 

sulating state, and the ADW state. To investigate these 
states, we calculate the superfluid order parameter $ = 
{ciaiCiai ), the quasiparticle weight Z, and the ADW or- 
der parameter AfADW = (1/4) (ns — nX), where n B ^ — 
J2 a ,a( n B(A) a a) with n B(A)aa = c B{A)mij c B{A)a<J being 
the number operator of the B(A) sublattice per site. Z 
represents the coherent spectral weight of the Bogoliubov 
quasiparticle (23[. Because of particle-hole symmetry, $ 
is independent of a. 

We first calculate Z and <& for the investigation of 
the superfluid-Mott insulator transition, without divid- 
ing the system into two sublattices. In Fig. QJa), the re- 
sults for D = are shown. As \U\ increases, Z decreases 
and jumps to with vanishing <I> simultaneously. The 
results indicate that the discontinuous quantum phase 
transition from the superfluid to the Mott insulator oc- 
curs. This Mott transition is caused by the two-band 
effects [lfj . We find that the Bogoliubov quasiparticle is 



dzp(z)G a (z,uj), 



UJ + fl- (-1)' 
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FIG. 1: (Color Online) (a) The quasiparticle weight Z and 
the superfluid order parameter $ as functions of U(< 0) for 
D = 0. (b) The ADW order parameter A/adw and Z as 
functions of U(< 0) for D = 0. Inset: The atomic number 
per site of the B(A) sublattice n B{A) = Y^ a , tT ( n B(A)c,a) ■ 



renormalized significantly towards the transition point. 

We next calculate Madw and Z, dividing the system 
into two sublattices. The results are shown in Fig. DJb). 
Z is nearly equal to 1 irrespective of U and Madw in- 
creases monotonously towards 1 with \U\. The results 
indicate that for D = the ADW order is enhanced by 
the attractive interaction U. As shown in the inset of Fig. 
[]Jb) the fillings of the neighboring sites approach 4 and 
0, respectively, with increasing \U\. The results demon- 
strate that the imbalanced atomic numbers between both 
sublattices take place in the ADW state and for \U\ > 4 
the almost fully-occupied and empty states emerge alter- 
nately. Note that ub + tia = 4 irrespective of U nor D 
due to particle-hole symmetry. 

To determine the most stable state among these three 
states, their energies are compared. As in Fig. H[a), we 
confirm that the ADW state is the most stable for all U 
at D = 0. For \U\ < 1 the energy of the ADW state 
is lower than that of the superfluid state, although their 
differences are so small in the present scale. 




FIG. 2: (Color Online) The energies of three states as func- 
tions of (a) U(< 0) for D = (b)D for U = -2.5. MI indicates 
the Mott insulator. 



We calculate Z, and Madw for D ^ by the same 
method. To investigate the most stable state among 
them, we compare the energies of these states. In Fig. 
[^b), the results for U = —2.5 are shown as functions 
of D. As D increases, the energy of the superfluid state 
decreases and crosses the energy of the ADW state. Ac- 
cordingly, the discontinuous quantum phase transition to 
the superfluid state occurs at D = 1.27. 




FIG. 3: (Color Online) The single-particle excitation spectra 
pAi(w) for several values of D at U = —2.5. 



We investigate the quantum phase transition from 
the spectral point of view. The single-particle ex- 
citation spectra (SPES) are defined as p ja (Lu) = 
— (l/Tr)ImG a (ui + i6) (7 = A, B), where S is a small pos- 
itive number, and psai^) and pAa{^) are given by the 
(1,1) and (2, 2) components of the matrix G a (u> + iS), re- 
spectively. Because of particle-hole symmetry, pai{oS) = 
PB2(—u) and pBi(to) = PA2(—^>)- For the superfluid and 
band-insulating states, the relation pA a {u) — pBa{u) is 
satisfied. In Fig. [3l the SPES pai(w) are shown for sev- 
eral values of D at U = -2.5. The SPES for D = 1.0 
shifts to the low energy region as compared to that for 
D = with scarcely changing the spectral shape. Ac- 
cordingly, the atomic numbers ua and ub for given U 
are independent of D. These properties are characteris- 
tic of the ADW state. The spectral gap around to = is 
~ 2.3, which is larger than D. The energy of the ADW 
state is thus independent of D for the range shown in 
Fig. [21^b). For D = 1.5, wc find the incoherent spectral 
weights away from u> = 0, which indicate the significant 
renormalization of the Bogoliubov quasiparticle. As D is 
increased, the enhanced incoherent spectral weights be- 
come inconspicuous and for D = 3.0, the SPES exhibits 
a typical profile of the band insulator. These findings are 
consistent with the behavior of other quantities. 
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We have confirmed for the single-band attractive Hub- 
bard model (J = U' = D = 0) at half-filling that the 
superfluid-insulator transition never occurs as shown in 
the studies so far [H, [H, [27[ . At half-filling the superfluid 
state and the ADW state are degenerate, while in the fill- 
ing deviated from the half-filling the superfluid state is 
the most stable [23|, [24], HEI . In the present system, the 
ADW state persists up to a certain value of D for given 
U. This feature is caused by the two-band effects. 



(a) 



(b) 




FIG. 4: (Color Online) (a) Quasiparticle weight Z as a func- 
tion of U(> 0) for D = 0. (b) The ADW order parame- 
ter Madw and Z at D = 0. Inset: The atomic number 
per site of the A sublattice with a (a) pseudospin riAa{a) = 
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FIG. 5: (Color Online) The energies of three states as func- 
tions of (a) U{> 0) for D = and (b)D for U = 2.0. MI 
indicates the Mott insulator. 



For the repulsive interaction region U > 0, we cal- 
culate Z and the staggered magnetization of the sub- 
= (1/2) (n^ - n lS ) (7 = A, B) with 



lattices -^^stagg — \ x l* i )\""y& 

n -ycr = XiaKau) (7 = A, B) and a being the opposite 
state of a. The relations M s ^ agg 



-M? 

stagg 



Maf and 

nA(B)o + nA(B)s = 2 are satisfied due to particle-hole 
symmetry. We first calculate Z for the investigation of 
the Mott transition, so that the system is not divided 
into two sublattices. As shown in Fig. 0Ja), the discon- 
tinuous transition from the metallic state to the Mott 
insulator occurs at U = 2.8. We next calculate Maf and 
Z, dividing the system into two sublattices. As shown 
in Fig. S£b), Maf increases toward 1 and Z takes the 
value close to 1 with increasing U. The inset of Fig. HJb) 



shows that the pseudospin state of the A sublattice ap- 
proaches the fully-polarized one, as \U\ is increased. The 
behavior exhibits typical features of the AF state. We 
compare the energies among the Mott insulating state, 
AF state, and renormalized metallic state to determine 
the most stable state. As shown in Fig. [SJ the AF state 
is the most stable for D = and the discontinuous phase 
transition to the metallic state occurs at D = 1.46 for 
U = 2.0. For U < 1 the energy of the AF state is lower 
than that of the metallic state, although their differences 
are so small in the present scale. 




FIG. 6: (Color Online) Phase diagram for the two-band Hub- 
bard model at half-filling. The spheres with blue (dark) color 
represent atoms with 'up' pseudospin and the spheres with 
red (light) color represent atoms with 'down' pseudospin. 



Systematically changing D and U, we investigate the 
ground state of the two-band Hubbard model at half- 
filling. The results are summarized in the phase diagram 
shown in Fig. [6] In the large D region, the band insula- 
tor comes into existence [10j . In the small D region, the 
ADW state and the AF state appear. As D is increased 
for the attractive region the direct transition between the 
ADW state and the band insulating state is observed for 
\U\ > 3.5, while for the repulsive region the metallic state 
appears between the AF and band insulating states even 
in U = 6. It is considered that the features are caused 
by the difference in spin and orbital fluctuations. In the 
repulsive interaction region, spin and orbital fluctuations 
are enhanced as U is increased [H, [2t| , leading to the sta- 
bilized metallic state even for large D. In the attractive 
interaction region, on the other hand, spin fluctuations 
disappear because of the formation of a Cooper pair. In 
the ADW state, orbital fluctuations also vanish. The 
ADW state is thus considered to be robust, which yields 
the direct transition to the band insulating state with 
increasing D. 



IV. DISCUSSIONS 

We investigate the stable ADW state using an effec- 
tive boson model, which is appropriate for the strong 
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attractive U in the half-filling. To derive the effective 
boson model from the Hamiltonian (jTJ) , we first make a 
particle-hole transformation in the higher orbital. We 
then expand the transformed Hamiltonian in the strong 
attractive U (D <C \U\), where two fermions with dif- 
ferent pseudospins make a hardcore boson. Within the 
fourth order expansion, the following effective model for 
the bosonic fermion pair can be derived, 

Woff = t cB &J,<Aa + U cff Y n ha n hl3 

(i,j),a i,a=£0 

+ Kff Y n i>°> n iSi ( 4 ) 

where bi <a annihilates a bosonic fermion pair on or- 
bital a in the ith lattice site, t cS (= —2t 2 /U) repre- 
sents the effective hopping integral of the boson, and 
Tii, a — b\ a bi, a . The effective interactions between two 
bosonic fermion pairs in the same lattice site and in the 
neighboring sites are denoted as U eS (= — 2U' + J — D) 
and V oS = — 2t^ s /U e g, respectively. Under the condi- 
tion U' = J = U/2(< 0), we obtain U cS > and thus 
V cS < 0. The repulsive U cS and attractive V eS prefer the 
ordered state where one bosonic fermion pair occupies 
each site with two orbitals being occupied alternately in 
the neighboring sites: the orbital AF state of the bo- 
son. Note that the occupied (unoccupied) state of the 
boson in the higher orbital represents the empty (doubly- 
occupied) state of the fermionic atoms in the original 
Hamiltonian, because of a particle-hole transformation. 
Accordingly, the orbital AF state of boson represents that 
the two orbitals in the same site are occupied by two 
pairs of fermionic atoms with different pseudospins and 
the neighboring site is empty in a viewpoint of the orig- 
inal Hamiltonian. This state is nothing but the ADW 
state. We have thus demonstrated that the effective bo- 
son model yields the ADW state in agreement with the 
numerical results. 

The effective Hamiltonian Q further suggests that for 
J = 2U' and D = the effective repulsion vanishes and 
Vjg is strongly enhanced. At these parameters, the su- 
perfluid state and/or the ADW state are expected to 
appear. To see the ordered state when the system ap- 
proaches the parameters J = 2U' and D — 0, we calcu- 
late the order parameters of the superfluid state <& and 
the ADW state Madw> releasing J from the restriction 
V = J = U/2. In Fig. [Jj the results for several values 
of J in U' = U/2 for D = are shown as functions of U. 
We find that $ keeps nonzero values only for J = 2U'. In 
this case, Madw/2 and $ exhibit the same U dependence 
as shown in Fig. [Jja). Furthermore, the energies of the 
superfluid and ADW states are degenerate at J = 2U' 
and D = 0, although the energy of the ADW state is 
the lowest for other parameters as shown in Fig. UJb). 
These results indicate that both states coexist: a super- 
solid state emerges. For the system with U' = J = U/2, 
the coexisting region of the superfluid and ADW states 
lies along the boundary of both phases shown in Fig. [S] 
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FIG. 7: (Color Online) (a) The ADW order parameter Madw 
and the superfluid order parameter $ as functions of the at- 
tractive interaction U for several values of J in U' = U /2 
and D = 0. For J = 2U' , Madw/2 and 4> show the same 
U dependence, (b) The energies of the ADW, superfluid, 
and Mott insulating states are compared for several values of 
J in U' — U/2 and D = 0. At the arrows, the superfluid 
(SF) states change into the Mott-insulating (MI) states. For 
J = 2U' , the energies of the superfluid and ADW states are 
degenerate. 



In real systems, there is a confinement potential. Re- 
cently, the effects of a harmonic confinement potential 
were investigated for the single-band attractive Hubbard 
model on the square lattice [l4[. It was shown that a 
harmonic confinement potential plays an essential role in 
stabilizing the supersolid state and that the doughnuts- 
like region of the supersolid state emerges between the 
ADW and superfluid phases. In the present system, the 
ADW and band-insulating states appear in addition to 
the superfluid state owing to the two-band effects. This is 
in contrast with the phases of the single-band attractive 
Hubbard model without a confinement potential, where 
only the superfluid state appears at half-filling as men- 
tioned before. For the two-band attractive lattice fermion 
systems with a confinement potential, a richer phase dia- 
gram including a supersolid state is expected to appear. 
The investigation of this issue with large numerical cal- 
culations is our next study. 



V. SUMMARY 

We have investigated the ordered state of ultracold 
fermionic atoms in three-dimensional optical lattices 
at half-filling using two-site dynamical mean-field the- 
ory. Because of the two-orbital effects, the ADW state 
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emerges for the attractive interaction region, while the 
AF state emerges for the repulsive interaction region. We 
have shown that the ADW and AF states are more sta- 
ble than the Mott insulating states. However, the ADW 
and AF states may become unstable by thermal fluctu- 
ations. Actually, it was shown for the single-band Hub- 
bard model that the transition from the AF state to the 
Mott insulator occurs around T/T-p ~ O(10 _1 ) at half- 
filling with increasing temperature (l6| . Here Tp is the 
Fermi temperature. The experiments for cold fermionic 
atoms in optical lattices have been performed in temper- 
atures T/T-p ~ C^IO" 1 ) %%. Accordingly, the compe- 
tition between the Mott insulating state and the ordered 
state (the ADW state or the AF state) at finite tempera- 
tures is expected. Detailed investigations for the ordered 
state of multiband fermionic systems at finite tempera- 
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